Weak lensing and dark energy: the impact of dark energy 
on nonlinear dark matter clustering 



Shahab Joudaki, Asantha Cooray 
Center for Cosmology, Dept. of Physics & Astronomy, University of California, Irvine, CA 92697 



o 

G 

in 

(N 



O 

U 

Oh- 

6 



> 

ON 
NO 
^1" 

o 

On 
O 



13 



Daniel E. Holz 

Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545 

(Dated: January 25, 2010) 

We examine the influence of percent-level dark energy corrections to the nonlinear matter power 
spectrum on constraints of the dark energy equation of state from future weak lensing probes. We 
explicitly show that a poor approximation (off by > 10%) to the nonlinear corrections causes a 
> la bias on the determination of the dark energy equation of state. Future weak lensing surveys 
must therefore incorporate dark energy modifications to the nonlinear matter power spectrum ac- 
curate to the percent-level, to avoid introducing significant bias in their measurements. For the 
WMAP5 cosmology, the more accurate power spectrum is more sensitive to dark energy properties, 
resulting in a factor of two improvement in dark energy equation of state constraints. We explore 
the complementary constraints on dark energy from future weak lensing and supernova surveys. 
A space-based, JDEM-like survey measures the equation of state in five independent redshift bins 
to ~ 10%, while this improves to ~ 5% for a wide-field ground-based survey like LSST. These 
constraints are contingent upon our ability to control weak lensing systematic uncertainties to the 
sub-percent level. 



I. INTRODUCTION 

The images of distant galaxies are gravitationally 
lensed by matter inhomogcncitics along the line-of- 
sight. In the weak lensing regime these percent- level 
magnifications and shape distortions of galaxies need 
to be analyzed statistically (see [HQ for a review). By 
extracting the shear power spectrum of weakly lensed 
sources [3h8| , the nature of the dark energy has been 
constrained with lensing surveys 0, Q ■ 

In a comprehensive analysis of future dark energy 
probes by the Dark Energy Task Force (DETF), weak 
lensing is singled out as particularly promising, in 
comparison with supernovae (SNe), galaxy cluster 
counting, and baryon acoustic oscillations (BAOs) @. 
An important aspect is that the lensing power spec- 
trum depends on both the lensing kernel and the 
matter power spectrum, making lensing a powerful 
probe of both background cosmology and the growth 
of structure. 

The optimism associated with lensing is predicated 
on overcoming the vast systematic uncertainties in 
both measurement and in theory [Iol - |22l |. These sys- 
tcmatics include dark energy corrections to the mod- 
cling of the nonlinear matter power spectrum [l8l . [l9j , 
higher order correction terms in the lensing integral 
(such as due to the Born approximation and lens- 
lens coupling [2(il - l22l | ) . and uncertainties of the mat- 
ter power spectrum on nonlinear scales due to the 
strong influence of baryonic physics p^4l7j . Ob- 
servational systematics include photometric redshift 
uncertainties, multiplicative factors in shear due to 
calibration errors, and additive factors due to PSF 
anisotropies (lOl. 

Furthermore, the observed cllipticities of weakly 
lensed galaxies are sensitive to the reduced shear, 



g = 7/(1 — k), where 7 is the shear and n is the 
convergence. In the weak lensing regime we make use 
of the expansion of the reduced shear to first order 
in the fields: g ~ 7. For future lensing surveys, this 
approximation induces a bias on the cosmological pa- 
rameters at the same order as that of the parameter 
constraints [23l, |24|. For current purposes, however, 
we continue to make use of this simplifying assump- 
tion of the shear as the lensing observable. 

In this work we examine one particular lensing sys- 
tematic: dark energy corrections to the nonlinear mat- 
ter power spectrum, and the impact of these on dark 
energy constraints from the weak lensing power spec- 
trum. We utilize two approaches towards modeling 
the nonlinear matter power spectrum in an evolv- 
ing equation of state (EOS) environment. Both ap- 
proaches are based upon the Smith et al. (2003) [25| 
prescription, valid for w = — 1. The conventional 
route in computing the matter power spectrum for 
w(z) =/= —I has been to use the w = — 1 fitting func- 
tions of Smith ct al. with appropriate modifications 
to the growth function and the redshift dependence of 
the matter density. We compare this method to that 
developed by McDonald, Trac, & Contaldi (2006) [H 
for a constant w cosmology, where numerical simula- 
tions underlie a fitting scheme that provides correc- 
tions to the Smith ct al. results. This latter ap- 
proach approximates the dependence of the matter 
power spectrum on a constant dark energy EOS to 
the level of a few percent. We further analyze mea- 
surements of the dark energy EOS in (decorrelated) 
redshift bins, as well as direct measurements of a two- 
parameter Taylor expansion form for the EOS. 

Our calculational methods are described in Sec- 
tion 2, with a basic review of the computation of 
distances, the growth function, the nonlinear matter 



power spectrum, and the weak lensing power spec- 
trum. In Section 3, we provide the dark energy con- 
straints from weak lensing tomography, utilizing two 
different approaches for calculating the nonlinear mat- 
ter power spectrum in a w{z) ^ — 1 cosmology. We ex- 
amine the complementarity between future weak lens- 
ing surveys and expansion history probes (e.g. super- 
novae distances measurements at z < 1.8). We also 
discuss the bias in dark energy due to uncertainties 
in the nonlinear matter power spectrum. At the end 
of Section 3 we provide an exploration of constraint 
contamination due to observational systematic uncer- 
tainties. Section 4 concludes with a discussion of our 
findings. We take our fiducial cosmological model, in 
accordance with WMAP5 data, to be a flat ACDM 
universe with fl c = 0.215, fl b = 0.045, h = 0.72, 
n s = 0.96, as — 0.8, and no massive neutrinos [26| . 



II. CALCULATIONAL METHOD 

Wc begin with a summary of our calculation. Wc 
briefly describe our cosmological distance and growth 
function, and then discuss the relevant observational 
quantities for weak lensing. 



Distances 



The comoving distance to an object at redshift z is 



1 



dz' 



(1) 



where X(x) = sin(a;) for a closed universe, sinh(a;) for 
an open universe, and x for a flat universe, and where 
Hq = 100 h km s _1 Mpc -1 is the Hubble constant. 
The Hubble parameter is given by 

h(z) = H oy /n m (i + z)3 + n w F( z ) + n k (i + z) 2 , 

(2) 

where {f2 m , £l w , f^} are the present matter, dark en- 
ergy, and curvature densities, in units of the crit- 
ical density. The evolution of the dark energy is 
represented by F(z). For a cosmological constant, 
F(z) = l. 

We parameterize the evolution of the dark en- 
ergy in two distinct ways. We use a popular Tay- 
lor expanded form for the dark energy EOS, given 
by w(z) = wo + z/(l + z)w a [27| , from which one 
obtains F(z) = (1 + z )^ 1+Wo+w ^e- 3Wa ^ . Alter- 
natively, instead of assuming a particular (physically 
unmotivated) model for the EOS, we utilize an agnos- 
tic, model-free approach to the redshift evolution of 
the dark energy [23 - l32j . We fit for w(z) binned in 
redshift: 



JY 



w(z) = ^2 w l E(z l ,z l+ i), 



(3) 



i=j=1 , w=-1 

- i=j=2, w=-1 
-■ - i=j=3, w=-1 

— i=j=4, w=-1 
i=j=5, w=-1 

- - i=j=5, w=-1/2 
■ ■ ■ i=j=5, w=-3/2 




FIG. 1. Geometric factor, Wi{z)Wj{z)x(z)/H(z), of five 
tomographic bins for our fiducial model (w = —1). For 
the fifth bin we also plot the kernel for w = —1/2 and 
w — —3/2. The narrowing of the lensing kernel for w > — 1 
stems from the decrease of each term with increasing w. 



where S(z^, Zi+i) is a tophat function in the region 
spanned by {zi, Zi+i}, but the ensuing analysis decor- 
relates the redshift bins. For this model-independent 
parameterization of the EOS, 



n-l 



FiZn-i <z<z n ) = (l+z)^ 1+w ^ Y[(l+z t ) 

i=0 



B. Growth Function 



3(wi— Wj+i) 



(4) 



For a particular parameterization of the EOS one 
can calculate the growth of matter fluctuations in the 
universe. For matter perturbations on linear scales, 
it is possible to separate out the time evolution of 
the perturbation: <5(k, z)/5(k, 0) = D(z)/ D(0), where 
D(z) is the growth function (which evolves as the 
scale factor in a matter dominated universe). For 
w(z) 7^—1 this linear approximation breaks down on 
very large scales due to clustering in the dark energy. 
However, due to the large uncertainties from cosmic 
variance at these scales, the impact of dark energy in- 
homogeneities are negligible for weak lensing studies 
(e.g. (SSI). 

The normalized growth function, G(z) = (1 + 
z)D(z), can be found by solving a second-order dif- 
ferential equation [34| : 
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This differential equation is valid for non-flat geome- 
tries, and carries the initial conditions {G(z m d) = 
1; ^| 2md = 0} in a matter dominated epoch z m( j. 
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C. Weak Lensing Observables 

We employ weak lensing tomography, wherein we 
divide the redshift distribution of source galaxies into 
distinct redshift bins 35 1. This provides informa- 
tion about the redshift distribution of the interven- 
ing lenses, and thereby allows for more stringent con- 
straints on cosmological parameters [HI HH . 

The number density of source galaxies in a square 
arcminute in each tomographic redshift bin (with 
boundaries < z s < Zi+i) is defined by = 
J^ +1 dz s p(z s ), where p(z s ) = ngj^e~^ z ' /zo ^ is the 

redshift distribution of source galaxies 3 . We adopt 
{zo — 0.5, a = 2, 13 = 1} [Ifj, appropriate for the 
Large Synoptic Survey Telescope (LSST [33,113), nor ~ 
malized such that J a dz p(z) = n g . We use the same 
distribution to describe our Joint Dark Energy Mis- 
sion ( JDEM) source population, but with improved n g 
(see Table QJ. We use the Born approximation, and 
perform the lensing calculation along the unperturbed 
photon path [20l - [22| | . The lensing weight function of 
i th tomographic bin is given by 



the " th 



W t (z) = 1^1(1 + z) X (z) x 



2 Tli 

21+1 dz A x{z :):^ z)) P {; 

m&x.(z.Zi) 



(6) 



for z < Zi+i, where x(z) = VWk\ Jo TT^jm ■ To 
ensure that we have no lenses behind our sources, we 
take Wi(z) = for z > z i+ \. The weight function 
increases for a more negative dark energy equation 
of state, as shown in Figure [1] The power spectrum 
of the convergence field is subsequently given by the 
Limber approximation [39j : 
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FIG. 2. Ratio of the present matter power spectrum of a 
constant dark energy EOS to that of a cosmological con- 
stant, A^ u (k, z)/A.\(k, z), as a function of w. The solid line 
is obtained by modifying the growth function and matter 
density evolution based on the Smith et al. prescription, 
while the other (more accurate) curves utilize the MTC 
corrections to the power spectrum. 



LSST [36|, |37j. We assume that the survey covers 
half of the sky, with a galaxy density of n g = 50 
30L 37 1 . For comparison, we also consider 



a space-based survey, such as a JDEM candidate like 
the SuperNova Acceleration Probe (SNAP) |43ti45|. 
For simplicity, we keep the same source distribution 
and intrinsic shear uncertainty, modifying the source 
density to twice that of the ground-based survey, and 
the width of the survey to a tenth of the sky. The 
characteristics of the two surveys are summarized in 
Table [J 
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W l (z)W J (z) X (z) 



dZ , A &l(*/x(*),*), 

(7) 

where zjj is the horizon redshift, and A^j L (Z/x(z), z) 
is the full nonlinear matter power spectrum. 

For I > 300, the dominant contribution to the mat- 
ter power spectrum will be from nonlinear scales [40l - 
HH, emphasizing the need to correctly model the ef- 
fect of dark energy at these scales. The observed con- 
vergence power spectrum, which is identical to that 
of the shear Q, is contaminated by shot noise due 
to the finite source density, as well as uncertainty 
in the intrinsic shapes of the source galaxies, lead- 
ing to: Cij(l) = Cij(l) + 5ij (j 2 ) jfii (which assumes 
that the noise is uncorrelated between tomographic 
bins). Whereas cosmic variance dominates the error 
on large angular scales, the shot noise is dominant on 
small scales. 

For simplicity we take the intrinsic shape uncer- 
tainty of the source galaxies to be redshift indepen- 

1 /2 

dent: ( 7 2 ) 1/z = 0.22, in accordance with expected 
results from a future ground-based survey such as 



D. Nonlinear Matter Power Spectrum 

We now detail our calculation of the matter power 
spectrum for a general dark energy cosmologyQ In the 
linear regime the transfer function T(fc, z) is computed 
using the prescription of Eisenstein & Hu (1997) [13]. 
The dimensionless linear power spectrum, normalized 
to the variance of the matter density field on scales of 
8 h^ 1 Mpc, er 8 , is given by 



A? 



(k,z) 



*T 2 (k,z) 



C8,obs 



D 2 (0) Wth 



(8) 



1 We do not pursue a direct halo-model approach |4(| for mod- 
eling the weak lensing power spectrum, since for the case of 
general dark energy models we have inadequate descriptions 
for the halo mass function, the halo dark matter profile, and 
the large-scale halo bias. 
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FIG. 3. Ratio of the matter power spectrum of a constant 
dark energy EOS to that of a cosmological constant, at 
z=0.75. 

where n s is the spectral index, <7s,obs 

= 0.8, and 



f 8, theory 



dlnfcA|(fc,0)J 2 (8fc') 



(9) 



with the spherical tophat filter, J(8k') = 
[3/(8fc') 3 ][sin(8fc') - {8k') cos(8fc')], and k' = k 
h^ 1 Mpc is dimcnsionless. One can subsequently 
extend this power spectrum to nonlinear scales by 
calculating the appropriate effective spectral index, 
effective spectral curvature, and nonlinear scale, 
em ploy ing the fitting functions provided in Smith et 
al. [If. 

The underlying cosmology in the Smith et al. fit- 
ting functions manifests itself in two distinct ways. 
First, cosmology impacts the evolution of the matter 
density, Q m {z), and the evolution of the growth of 
matter perturbations, D{z). Second, cosmology fixes 
the functional form and coefficients associated with 
the fitting functions, which are fine-tuned to a suite 
of ACDM N-body simulations. Thus, whereas an ar- 
bitrary dark energy EOS could make an imprint on 
the matter power spectrum via its influence on £l m {z) 
and D{z), the cosmological dependence of the N-body 
fitting functions remain fine-tuned to a w = — 1 dark 
energy EOS. A correct approximation of the matter 
power spectrum in an evolving EOS cosmology there- 
fore requires both an update of the effective spectral 
parameters of the Smith et al. prescription, and a gen- 
eralization of the fitting functions. 



Probe 


/sky 


n g (arcmin 2 ) 


■^pcak 




LSST 


0.5 


50 


1.0 


0.22 


JDEM 


0.1 


100 


1.0 


0.22 



TABLE I. Descriptions of our fiducial ground-based 
(LSST) and space-based (JDEM) probes. 



No. of bins 


2l Zo Z'i Za Z<\ 


1 bin 


3.0 - 


2 bins 


1.3 3.0 - 


3 bins 


1.0 1.6 3.0 


4 bins 


0.85 1.3 1.85 3.0 - 


5 bins 


0.75 1.1 1.45 1.95 3.0 



TABLE II. Redshift bins of the source distribution. The 
redshifts are determined such that each tomographic bin 
contains roughly the same number density of galaxies. We 
ignore sources above z = 3. 



A full calculation of the matter power spectrum will 
thus require a large suite of N-body simulations. An 
alternative approach devised by McDonald, Trac, & 
Contaldi (2006; hereafter MTC) notes that many nu- 
merical uncertainties cancel when taking the ratios 
of power spectra [lj . By interpolating the matter 
power spectrum ratios between constant EOS values, 
MTC provide a route to a fast calculation of the non- 
linear matter power spectrum in constant w cosmolo- 
gies. An accurate matter power spectrum can thereby 
be obtained by computing the power spectrum for a 
ACDM universe, and subsequently multiplying it by a 
cosmology-dependent correction factor. In particular, 



D 2 {z,w)/D 2 {0,w) 
D 2 {z,w)/D 2 {0,w) 



Anl,a(M) x 



(10) 



T(fc,z,p) 



where the N cosmological parameters are given by 
p = {Q m ,ttb,h : as,n s ,w}, w = -1, and 



T{k,z,p) 



,(fc) 



(11) 



In this equation, M is the polynomial order of the 
parameters, and the coefficients A VlU2 ... UN {k) are pro- 
vided by a least squares fit to the MTC simulations. 
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FIG. 4. Ratio of the matter power spectrum of a constant 
dark energy EOS to that of a cosmological constant, at 

55=1.5. 
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FIG. 5. Convergence power spectra for the case of five 
tomographic bins in a ACDM cosmology. We also include 
a number of representative cross-terms. For the i = j = 5 
case, we include the noise contribution as a band about 
the curve. 

On large scales the correction factor is equal to unity, 
as expected. The power spectrum ratios are only pro- 
vided up to z = 3/2 and k = 10 h Mpc -1 . How- 
ever, since we only consider multipoles up to I = 3000, 
scales smaller than k — 10 h Mpc -1 are only probed 
at very low rcdshift, where the number density of 
sources becomes negligible. Moreover, as gradually 
smaller scales become linear with increasing redshift, 
for a fixed range of scales, the influence of dark energy 
on the matter power spectrum is progressively mod- 
elled by linear theory as the redshift increases. The 
dark energy correction factor for the matter power 
spectrum therefore covers an adequate range of scales 
and rcdshifts for weak lensing studies. Furthermore, 
as the MTC simulations have only been carried out 
for a limited range of parameter space, we have en- 
sured that all of the cosmological parameters under 
consideration here live within their explored region. 

III. RESULTS 

We next explore the differences between the Smith 
et al. and MTC matter power spectra, and their re- 
spective impacts on cosmological constraints from fu- 
ture weak lensing probes. Systematic uncertainties 
are discussed in Section G. 



A. Matter Power Spectrum: MTC vs. Smith 

In Figures 2-4 we show the evolution of the MTC 
nonlinear matter power spectrum as a function of dark 
energy EOS (for constant w). Each figure represents 
a different rcdshift (z = 0, 0.75, and 1.5), and every 
plot shows three different MTC curves, corresponding 



to scales of k = [0.1, 1.0, 10] h Mpc -1 . The curves are 
normalized to the power spectrum corresponding to a 
pure cosmological constant, and the Smith et al. fit 
(modifying solely the growth function and the matter 
density in accordance with the dark energy EOS) is 
also shown, for comparison. The MTC power spec- 
tra show large (> 10%) deviations from the simplified 
Smith et al. approximation over a range of reasonable 
w values. In general, constant w < — 1 models lead to 
a suppression of the matter power spectrum as com- 
pared to that in a ACDM universe, with the suppre- 
sion more pronounced at smaller (nonlinear) scales. 
Conversely, the matter power spectrum is enhanced 
for w > — 1. As expected, the MTC correction factor 
becomes less pronounced at higher redshift as increas- 
ingly smaller scales become linear. 

B. Weak Lensing 

The convergence power spectrum depends on both 
the geometric factor, Wi(z)Wj(z)x(z)/ H(z), and the 
matter power spectrum, A 2 (k,z), but in a compet- 
ing manner. Whereas a decreasing dark energy EOS 
has the effect of increasing the geometric factor, and 
therefore the convergence power spectrum, it also sup- 
presses the nonlinear matter power spectrum, which 
in turn suppresses the convergence power spectrum. 
As a result, these two effects partially cancel, thereby 
decreasing the sensitivity of weak lensing to dark en- 
ergy. 

In Figure [5] we plot the power spectrum of the con- 
vergence for five tomographic bins (with the redshift 
divisions listed in Table HI]) , as well as some repre- 
sentative convergence power spectra cross-terms (e.g., 
between bins 4 and 5). We divide the power spectrum 

by its noise, ACy-(Z) = 2 ^2/ (21 + 1){C VJ (I) + 




FIG. 6. Signal to noise with five tomographic bins in a 
ACDM cosmology, for LSST (denoted by 'L') and JDEM 
('J'). 
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FIG. 7. Ratio of MTC to Smith et al. convergence spec- 
tra derivatives with respect to w, i.e. DMTc/Dsmith = 
dc, I /8Cn 
~l MTC /"e^rlsmith. 

8ij ^7 2 ) /fii), in Figure HJ The signal to noise is con- 
sistently higher for the wider LSST than it is for the 
deeper JDEM. The larger width of LSST gives strong 
signal to noise in particular at smaller multipoles, 
while the greater depth of a JDEM-type survey makes 
it increasingly competitive at larger multipoles. 

Figure [7] presents the ratio of the convergence spec- 
tra derivatives with respect to a constant EOS be- 
tween Smith et al. and MTC. The MTC derivatives 
approach the Smith et al. derivatives in higher-redshift 
tomographic bins, as the corresponding matter power 
spectra converge at large redshift. This is because the 
MTC correction factor approaches unity at high z (see 
Figs. 2-4), and we therefore impose the approximation 
A 2 (M)|mtc = A 2 (fc,z)| Srai th for z > 1.5. TablesEI] 
and[V]show that our conclusions are insensitive to this 
approximation. The z < 1.5 entries in these tables 
are for five tomographic bins, with redshift divisions 
at [0.59,0.83,1.04,1.26,1.5], engineered to preserve an 
equal number density of sources across each bin. Since 
the MTC derivatives have decreased relative to the 
Smith et al. ones, one might conclude that the sen- 
sitivity to dark energy has also decreased. However, 
due to the particular nature of the cross-correlations 
between cosmological parameters, this turns out not 
to be the case for our fiducial cosmology, as is shown 
below. 



C. Constraints on Dark Energy 

In the previous subsections we have examined the 
corrections to the nonlinear matter power spectrum 
due to the presence of dark energy. We now explore 
how these corrections impact weak lensing constraints 
of the dark energy. Assuming Gaussianity in the like- 
lihood distribution of the cosmologic al p arameters for 
the lensing power spectrum (l3l . [501 - 4531 ] , we utilize a 



Fisher matrix analysis [48], |4t| : 

V (2Z + 1)A/ [dCt-.dd-,- 

(12) 

We consider constraints for both a ground-based sur- 
vey, such as LSST, and a space-based JDEM probe, 
such as SNAP. We analyze multipoles between 50 and 
3000, with the power spectrum of the convergence 
field divided into multipole bins of width Al = 50, 
(the results are insensitive to the choice of binning). 
The cutoff at I = 3000 avoids non-Gaussianities of the 
convergence field [HI, I50T453T ] , as well as uncertainties 
from baryonic physics that increase at larger multi- 
poles mGi. 

We present constraints on cosmology for up to five 
tomographic bins, where each bin is constructed to 
contain the same effective number density of galax- 
ies. The redshift divisions of the bins are listed in Ta- 
ble [Til For the terms in Eq. [T^Jwe carry out two-sided 
numerical derivatives with steps of 10% in parameter 
value (except for ±0.1 in the case of w a , and ±0.05 
for fifc). Our results are essentially unaffected by the 
particular choice of step size. We note that there is 
an ambiguity in the evaluation of the flk derivative 
at Qk = 0, since the Smith et al. fitting functions are 
only given for flat and open universes. Although it is 
common to carry out the curvature derivative using 
the flat cosmology fitting functions, we also evaluate 
the derivative at a slightly open universe to see if it 
renders noticeable differences in the constraints. This 
assumes that the power spectrum doesn't vary sharply 
about flatness. The two methods for calculating the 
ilfc derivative yield consistent results (within ~ 15%), 
and we quote results utilizing the latter method. 

We consider two dark energy models: w(z) = wo + 
z/(l + z)w a , and w(z) = w c is constant (i.e., w a = 0). 
For the w c case, we define the figure of merit (FOM) 
of the EOS measurement to be FOM(w c ) = <r~ 2 (w c ). 
For the {uio,w a } case, we take the FOM to be the 
inverse of the error ellipse in the {wo,w a } plane (as 
defined in the DETF report @): 

FOMdetf = [a(w p )(j{w a )]-\ (13) 

where w p is the pivot value of the EOS. We do not in- 
clude the factor of 1/tt in the DETF FOM Note 
that the uncertainty on the dark energy EOS at the 
pivot is equal to the uncertainty on a constant EOS § . 
We find that there is roughly a factor of two improve- 
ment in FOM(w c ) in going from Smith et al. to the 
MTC form for the dark matter power spectrum, as 
shown in Table UTT1 This relative improvement in the 
EOS constraint is independent of the number of tomo- 
graphic bins. Although the derivatives in Eq. [12] are 
larger for the Smith et al. case (see Fig. [7]), the MTC 
form for the power spectrum provides more stringent 
constraints due to its reduced cross-correlations in the 
Fisher matrix. For a 2 x 2 matrix, this can be visu- 
alized as rendering a larger determinant, and thereby 
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0.021 
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140 


0.025 


1600 


2 (M-JDEM) 


0.22 


0.71 


31 


0.046 


480 


5 (M-JDEM) 


0.11 


0.44 


100 


0.023 
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+n k 


0.22 


0.96 


12 


0.025 


1700 


+ Priors 


0.11 


0.44 


98 


0.023 
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0.30 
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22 


0.036 


760 


+Q k 
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+ Priors 
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15 
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+ Priors 


0.015 


0.56 


57 
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1000 


+ 2300 SNe 


0.048 


0.21 


230 


0.021 


2200 


2 (S-JDEM) 


0.32 


0.79 


21 


0.059 


280 


5 (S-JDEM) 


0.11 


0.34 


87 


0.034 


860 


2 < 1.5 


0.14 


0.44 


51 


0.041 


580 



TABLE III. FOM values and la constraints on the dark 
energy parameters, based on either the Smith et al. pre- 
scription (denoted by 'S') or with the inclusion of MTC 
corrections (denoted by 'M'). The +2300 entries in the 
table are a joint analysis of WL and SNe, with SN sys- 
tematics included (flatness assumed). The 2300 SNe alone 
constrain [wo,w a ,Wc] to [0.04, 0.67, 0.036] in the case of 
statistical noise only, and to [0.055, 1.0, 0.054] when sys- 
tematic errors are included. For the cases with priors, we 
utilize an HST prior of 0.08 on h [55| and a Planck prior 
of 0.0032 on Q, k [56| . The strongest improvement comes 
from the curvature prior. 



better constraints as the inverse of the Fisher matrix 
is taken. Thus, even though smaller derivatives com- 
monly provide poorer constraints, here we find that 
the cross-correlation terms yield improved constraints 
for a constant (as well as a redshift-binncd) EOS. 

Generalizing to non-flat cosmologies significantly 
degrades the cosmological constraints, especially if we 
only consider effects out to z = 1.5. This degrada- 
tion is ameliorated with an expected Planck prior of 



0.0032 on fl k [56J and HST prior of 0.08 on h [55|. 
For LSST, using the Smith et al. prescription beyond 
z = 1.5, the inclusion of curvature drives the FOM 
down by one sixth from 4600 to 3800, whereas in the 
case of a redshift cutoff (at z = 1.5) the inclusion of 
the curvature parameter causes a factor of five deteri- 
oration of the dark energy FOM from 2900 to 550. By 
including HST and Planck priors, the FOM increases 
to 4500 and 2700, respectively. A similar sensitivity 
to the curvature density is seen for the JDEM case, 
which is also effectively removed by our choice of pri- 
ors. We moreover consider a lower cutoff in / at 1000, 
noting that this cutoff alone would diminish the con- 
straints by roughly a factor of two for LSST and 2.5 for 
JDEM. The relative constraint degradation is larger 
for a JDEM probe as its integrated signal-to-noise is 
more sensitive to larger multipoles (as seen in Fig. [6]). 



D. Evolving Dark Energy EOS 

One major limitation of the MTC correction to the 
nonlinear dark matter power spectrum is that it is 
only valid in the case of a constant dark energy EOS. 
We now generalize this approach to the case of dy- 
namical dark energy. We utilize an observation by 
Francis et al. [54[, who have shown that the power 
spectrum of a dynamical w(z) resembles that of a 
constant w with the same distance to the last scat- 
tering surface. In other words, when the integrated 
expansion history between a dynamical EOS and a 
constant EOS are similar, the growth histories will 
likewise be similar. Francis et al. have numerically 
confirmed this to a few percent, at scales up to k = 5 
h Mpc -1 for the {wo,w a } parameterization. Based 
upon this, we mimic the matter power spectrum of an 
evolving dark energy EOS by utilizing a constant EOS 
with a matched distance to the last scattering surface. 
We calculate all other components of the convergence 
spectrum based upon the actual w{z). 

The FOM of the DETF remains relatively un- 
changed despite the inclusion of the improved (MTC) 
nonlinear matter power spectrum. For the {wo,w a } 
parameterization, we note that although the con- 
straints due to changes in wq are improved for MTC, 
the Smith et al. constraints are more sensitive to 
changes in w a , and therefore their net effects on the 
FOM are roughly equivalent. This latter decrease in 
sensitivity for MTC occurs as the distance-matching 
prescription of mapping {wq = — l,w a } onto w c ren- 
ders larger deviations in the MTC matter power spec- 
trum, and thereby smaller deviations in the conver- 
gence spectrum due to the cancellation with the lens- 
ing kernel discussed above. The smaller w a derivatives 
then lead to inferior constraints. Imposing a cutoff at 
I = 1000 renders a factor of five deterioration in the 
FOM of the DETF, primarily due to the factor of 
three deterioration of the constraint on w a , as shown 
in Table HIU 
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In addition to the two-parameter {wo,w a } con- 
straints discussed above, we also consider constraints 
on a redshift-binned EOS. We utilize the MTC 
method, generalized following the Francis et al. ap- 
proach of matching the last-scattering distances be- 
tween evolving and non-evolving dark energy. Once 
the weak lensing Fisher matrix for correlated dark en- 
ergy bins is obtained, we rotate the dark energy pa- 
rameters into a basis where they are uncorrelated. 

The most straightforward approach would be to 
diagonalize the marginalized Fisher matrix F (i.e. 
marginalized over all parameters except for the EOS, 
Wi=i...N), such that the uncertainties in the uncor- 
related dark energy parameters are given by the in- 
verse of the eigenvalues. However, we choose the 
transformation matrix advocated in Huterer & Cooray 
(2005) dj, namely F 1 / 2 = T A 1 / 2 0, where O 
and A are the eigenvector and eigenvalue matrices. 
The uncorrelated parameters are then given by q = 
F 1/2 P [H, HI HI EH- This transformation matrix 
has the pleasing feature of generating relatively local- 
ized and mostly positive weights |3CJ, E3 , while also 
preserving the information content of the correlated 
Fisher matrix. Note that the transformation matrix 
is normalized such that the sum of the elements along 
each row is equal to unity [30|, EH . The rows therefore 
represent the contribution of each redshift bin to the 
decorrelated parameters. The errors on the uncorre- 
lated parameters are then given by 



(A qi A qj ) 



a p 



(14) 



We now define the figure of merit for this param- 
eterization to be that of Sullivan, Cooray, & Holz 
(2007) [p]: 



FOM 



■wbii 



''(Wi) 



W (15) 



We note that this FOM is equal to that of a constant 
EOS [59| . which is an important feature as it explic- 
itly demonstrates that the FOM is binning indepen- 
dent. This can be seen by comparing the case of seven 
bins, tUj = i...7, with that of one bin, w c , in Tables ITPTl 
and[V] In addition, as the parameter constraints scale 

— 1/2 

as / sky , both FOMwbinncd and FOMdetf increase 
linearly with the fraction of the sky covered. 

For comparison with the EOS constraints in Sarkar 
et al. , we have chosen the same redshift binning of the 
EOS (boundaries at z = {0.07, 0.15, 0.3, 0.6, 1.2, 3.0}). 
Although the seventh bin is highly unconstrained, the 
uncertainty does not leak noticeably into the other 
parameters (as compared to fixing = —1). This 
particular binning attempted to maximize the num- 
ber of redshift bins with 10% or better constraints on 
the dark energy EOS from future cosmic microwave 
background (CMB), BAO, and SN data. We have not 
provided a comparable binning for the case of weak 
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TABLE IV. 1<t uncertainties on cosmological parameters 
from WL alone, for the case of five tomographic bins in a 
flat universe without external priors. 'S' stands for Smith 
et al., and 'M' for MTC. 



lensing. One possible approach is to pick redshifts 
such that the sensitivity of the integrand to changes 
in w arc equivalent in each bin. 

Table IIVI shows the uncertainties on cosmological 
parameters for three different EOS parametcrizations 
(w c , {wo,w a }, and redshift-binned w(z)), as deter- 
mined from weak lensing with five tomographic bins. 
For a redshift-binned w(z), including the MTC cor- 
rections to the dark matter power spectrum, an LSST 
weak lensing survey constrains five EOS parameters 
to better than 10% (three of which are better than 
5%), and a JDEM-like survey constrains three EOS 
parameters. The LSST binned FOM is a factor of 2.5 
larger than the one for JDEM, and the constraints on 
other cosmological parameters are roughly a factor of 
two better. The main point here is that a very wide 
survey is more effective than a deep, but narrow, sur- 
vey. Although a deep, moderately wide survey could 
be even more effective, the constraints depend sensi- 
tively on the precise nature of the surveys. The MTC 
improvement to the nonlinear dark matter power spec- 
trum leads to a factor of two improvement in the FOM 
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No. of bins 


cr(ioi) fr(«)2) ^(^"3) ""(104) (7(105) ""(tog) FOM 


1 (M-LSST) 


0.59 1.47 0.96 0.69 0.74 1.6 8.8 


z ^IVl-ijoo 1 ) 
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0.083 0.074 0.037 0.020 0.032 0.33 4600 


I < 1000 
+ Priors 


1 \ 1 it r \ \ 1 \ n n f n a non n n r - n n 1 n onn 

0.12 0.10 0.052 0.029 0.050 0.16 2200 

n 1 n 1 ~i n nc n n notr n n^c n on 1 ,1 nn 

(J. lo U.ll U.UbU U.Uoo U.U75 U.zU 14UU 
niQ n 1 n n nco n hqq n n/i n. nin oofin 

u. io u. iu u.uoz u.uzo u.u4y u.iy zzuu 


z < 1.5 
in 

-f- Priors 


nil nin n n/i c; nno/i n n/i o ~\ a onnn 
U.ll U.1U U.U40 U.Uz4 U.U4o 1.4 oUUU 

n on n i a. n nac n n^o n /in n f^o f^fin 

U.ZU U.lo U.Uoo U.UOZ U.4U U.DZ ODU 

n 19 nin n nzLfi n 09^ n oar 1 97 97nn 
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0.055 0.049 0.032 0.021 0.032 0.50 5000 
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3 (S-LSST) 


0.33 0.40 0.12 0.037 0.035 0.43 1700 


4 (S-LSST) 


0.22 0.29 0.095 0.031 0.034 0.57 2100 


5 (S-LbS 1 ) 
z < 1.5 


n rr n o-i n nnn n n nn n n ,1 n tro oonn 

0.25 0.31 0.090 0.029 0.034 0.53 2200 
0.26 0.28 0.087 0.032 0.047 1.0 1600 


2 (M-JDEM) 


0.22 0.22 0.12 0.067 0.082 0.52 490 


5 (M-JDEM) 

+n k 

+ Priors 


0.12 0.11 0.056 0.031 0.050 0.52 1900 
0.12 0.11 0.056 0.032 0.062 0.22 1700 
0.12 0.11 0.056 0.031 0.051 0.49 1900 


£ < 1000 

+n k 

T Priors 


non /I it n no r* n n a a nnop" no/ - * '■znn 

0.20 0.17 0.08b 0.048 0.085 0.26 790 

n oo n i o n nn t n ncr nio nin rrtrn 

0.22 0.18 0.097 0.05b 0.12 0.32 550 
0.21 0.1/ 0.08b 0.04/ 0.081 0.30 810 


z < 1.5 
+ Priors 


0.18 0.16 0.073 0.040 0.078 1.7 1100 
27 26 13 075 67 75 270 
0.18 0.16 0.073 0.040 0.079 1.6 1000 


+ 2300 SNe 


0.067 0.061 0.045 0.033 0.056 0.58 2200 


2 (S-JDEM) 


0.45 1.3 0.43 0.11 0.075 0.29 280 


5 (S-JDEM) 
z < 1.5 


0.36 0.48 0.14 0.046 0.055 1.5 870 
0.36 0.45 0.14 0.052 0.082 1.5 580 



TABLE V. la constraints on values of w(z) in uncor- 
related redshift bins, where 'S' stands for Smith et al., 
and 'M' for MTC. The EOS is redshift binned at z = 
{0.07,0.15,0.3,0.6,1.2,3.0}, but decorrelating the covari- 
ance matrix for the original bins results in a leakage across 
bins. Figure [8] shows this leakage in terms of window func- 
tions. The 2300 SNe alone constrain the six EOS bins to 
[0.11,0.10,0.12,0.69,0.24,1.2]. 



of a rcdshift-binncd EOS. 

It is possible to compare our weak lensing con- 
straints with the combined constraints from SNe, 
BAOs, and CMB. The latter have been calculated in 
Sarkar et al. [32|, using a Markov Chain Monte Carlo 
(MCMC) likelihood approach. Case C in Sarkar et al. 
considers a mock catalog of 300 SNe at z < 0.1 
and 2000 SNe in the range of 0.1 < z < 1.8 0, two 
current BAO distance estimates [gj [62|, fifteen op- 
timistic future BAO measurements (five from V1N1 
of [H], ten from ADEPT (H), as well as CMB con- 
straints on [Q m h, h, R] of [0.023, 0.08, 0.03] [H, M, 



|66( , where R is the distance to the last scattering sur- 
face. For this extensive data set, the redshift-binned 
dark energy FOM is 1300. In comparison, our pro- 
jected lensing constraints yield an FOM of 4700 for 
LSST and 1900 for JDEM (for the MTC matter power 
spectrum). It is evident that weak lensing is poten- 
tially an extraordinarily powerful dark energy probe. 



E. SN/Weak Lensing Complementarity 

Our distribution of lensing source galaxies peaks 
at a redshift of unity. It is therefore interesting 
to combine our highcr-redshift lensing measurements 
with those from supernovac. To this end, we uni- 
formly distribute a set of 300 SNe at z < 0.1 [<30[, 
and 2000 SNe in the range 0.1 < z < 1.8 (as ex- 
pected from a space-based JDEM probe [Hj|, or as 
part of a first data release by a ground-based tele- 
scope such as LSST). For each supernova we take the 
intrinsic noise to be a Gaussian in magnitude with 
Cint = 0.1 [67|- We divide the Hubble diagram for 
z > 0.1 into 50 redshift bins, and associate each bin 
with a redshift-dependent systematic floor of magni- 
tude S m = 0.02(0. l/Az) 1 / 2 (1.7/z max )(l + z)/2.7 @. 
We assume that this irreducible systematic has no cor- 
relation between bins. For the SN constraints we use 
an HST prior of 0.08 for h [5|| (which is also applied 
to the WL case in the comparison with SNe+WL). 
For further details on the SN approach we refer the 




0.2 0.4 0.6 0.8 1.0 1.2 



FIG. 8. Lower portion of the figure shows the constraint 
on the dark energy EOS in five redshift bins. The narrow 
(black) error bars are obtained from weak lensing alone, for 
five tomographic bins, with an HST prior on the Hubble 
constant. The wide (red) error bars are obtained from a 
joint analysis of WL and SNe. Note that we are neglecting 
WL systematic uncertainties in this plot. The upper half 
of the figure illustrates the window functions associated 
with the first five decorrelated dark energy bins (WL + 
SNe). (Note that the window functions have been shifted 
down by a constant of 0.8 for visual clarity.) 
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TABLE VI. Bias on the determination of the dark energy 
equation of state, for the case of five tomographic bins (in 
a flat universe with no external priors). For the FOB of 
{wo,w a }, a value of 1.5 indicates a la shift from the true 
estimate, and a value of 2.5 indicates a 2a shift [23l. |70||. 
For w c , these values are 1.0 and 2.0, respectively. The 
lower entries list the bias on the correlated w(z) bins. 
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reader to earlier work [32 1 . 

The uncorrelated redshift-bin uncertainties in five 
tomographic weak lensing bins are given in Tables IIV 
and [V] and illustrated in Figure |5] for an LSST data 
set. The narrow (black) error bars are the constraints 
from weak lensing alone, while the wide (red) error 
bars represent the constraints from both weak lens- 
ing and SNe. Although weak lensing alone constrains 
three redshift bins to 5%, by combining this data set 
with SNe we constrain two additional EOS parameters 
to that level. Weak lensing and SNe are complemen- 
tary, as weak lensing is most effective for constraints at 
intermediate redshifts (roughly 0.2 < z < 1.2), while 
future SN data provide a more effective probe at lower 
redshifts (z < 0.2). 



F. Bias in dark energy due to uncertainties in 
the nonlinear matter power spectrum 

We have calculated the improvement in weak lens- 
ing constraints from use of the MTC matter power 
spectrum (incorporating w(z) ^= —1)- It is also inter- 
esting to calculate the bias that would arise from use 
of the more approximate Smith et al. power spectrum, 
when the real data is described by MTC. The bias in 
each parameter is given by [23l . Ir33 | : 



_ 1 (2Z + 1)A/ 



Tr 



Op/3 



5p a = /sky 

(16) 

where SCi is the difference between the MTC and 
Smith ct al. convergence spectra. One can calculate 
the corresponding figure of bias (FOB) for the subset 
of EOS parameters as rendered by [23| : 



1/2 



F0B = ( J2 s P a F a p6 Pfi 

,a,(8 



(17) 



Thus, for the case of one parameter, the FOB is simply 
equal to the ratio of that parameter's bias to its un- 
certainty. Table EJ presents the FOB values for LSST 



FIG. 9. Deterioration of constraints as a function of 
a Gaussian prior on each centroid of the redshift bins 
for LSST. We plot both FOM syst+ p r io r /FOM no _ syst and 
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FIG. 10. Deterioration of constraints as a function of the 
prior on the centroid of the redshift bins for a JDEM sur- 
vey. 



and JDEM weak lensing measurements. We find a la- 
2a bias between the Smith et al. and MTC determina- 
tions of the dark energy EOS (for a w = —0.9 fiducial 
cosmology). It is therefore critical that the nonlinear 
matter power spectrum be well-characterized, to en- 
able precision constraints on a dynamical dark energy 
equation of state. 



G. Systematic Uncertainties 

To assess the importance of systematic errors in the 
weak lensing measurements, we approximate a pos- 
sible shear miscalibration by including a multiplica- 
tive factor, fi, for each tomographic bin. We ignore 
additive systematic terms to the convergence power 
spectrum, as their dependence on I is currently poorly 
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FIG. 11. Deterioration of constraints as a function of the 
Gaussian prior on each multiplicative factor in shear for 
LSST. 



i — i — i i i 1 1 1 

JDEM 



_L_ 




-FOM Binned 



— w t 
■ — w, 

FOM DETF 

— Wo 
— 'W, 
- Ws 



1.0001 0.0010 0.0100 

Multiplicative Systematics Prior 



0.1000 



FIG. 12. Deterioration of constraints as a function of the 
prior on multiplicative factors in shear for a JDEM survey. 



understood [Til] - The uncertainty in the redshift dis- 
tribution of the sources is modeled to first order as a 
shift in the centroid of the tomographic bins [ll| . We 
ignore other possible sources of systematics, such as 
baryonic uncertainties in the nonlinear scales of the 

and intrinsic ellipticity 



14 



matter power spectrum 
correlations of the source galaxies 113. 

Following Hutcrer ct al. (2006) [ill ], the cross- 
correlation terms are changed by the systematics as: 
Cij (I) = Cij (I; z t + Szi, zj + Szj) [l + fi + fj}. For 
our five tomographic bins we therefore introduce an 
additional 5 x 2 = 10 free parameters to encapsulate 
possible systematics. The effect of these systematic 
uncertainties on the measurement of dark energy is 
heavily influenced by the corresponding priors. If we 
are optimistic, and claim an understanding of the sys- 
tematics at the 0.1% level, the deterioration in dark 
energy constraints is essentially negligible. However, 
the constraints are strongly compromised if the sys- 



tematics are unknown at the level of 1%. This is il- 
lustrated in Figures 9-12. 

Both of the surveys offer the potential for self- 
calibration tests, if the priors are above the percent 
level. For LSST, there is a factor of five (three) degra- 
dation of the binned FOM (DETF FOM) for multi- 
plicative shear errors and a factor of 14 deterioration 
in the two FOMs for errors in the redshift bin cen- 
troids. In general the degradations are slightly milder 
for JDEM; e.g. a factor of 11 degradation in the two 
FOMs for redshift bin centroid uncertainties. We note 
that a combined analysis with the lensing bis pec trum 
has the potential to improve these prospects [Tl| . 



IV. CONCLUSIONS 

We calculate weak lensing constraints on the dark 
energy equation of state, incorporating an improved 
nonlinear matter power spectrum which accounts for 
the effects of a dynamical dark energy. The most 
commonly utilized nonlinear extension of the matter 
power spectrum, calibrated from N-body simulations 
by Smith ct al., does not incorporate the effects of 
time-evolving dark energy. We follow a prescription 
presented by McDonald, Trac, & Contaldi (2006), and 
extend the Smith et al. form to include evolving dark 
energy. We then perform a full Fisher matrix analysis 
for two prospective weak lensing surveys, utilizing the 
improved power spectrum. 

By considering weak lensing tomography with an 
improved nonlinear matter power spectrum that in- 
corporates dynamical dark energy, we find a factor of 
two improvement in the dark energy figure of merit 
(for the WMAP5 fiducial cosmology). Although the 
changes in the matter power spectrum and the lens- 
ing kernel somewhat cancel in their effects on the con- 
vergence power spectrum, the parameter correlations 
nonetheless lead to improved parameter constraints. 

We further show that a poor approximation (off by 
> 10%) to the dark energy corrections of the non- 
linear matter power spectrum leads to a > ler bias 
on the dark energy equation of state. Future weak 
lensing surveys must therefore incorporate percent- 
level accurate dark energy modifications to the non- 
linear matter power spectrum to avoid introducing 
significant bias in their measurements. In addition, 
by combining weak lensing data with supernova mea- 
surements at lower redshifts, we show that a general 
dark energy model can be constrained in five redshift 
bins to 5% for an LSST-type survey and to 10% for 
a JDEM-like probe. These dark energy constraints 
are contingent upon our ability to understand weak 
lensing systematic uncertainties, such as those arising 
from shear miscalibration and the redshift uncertainty 
of the sources. If it is possible to control these sys- 
tematics to 0.1%, weak lensing constraints of dynam- 
ical dark energy from next-generation surveys offer 
tremendous promise. 
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